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Abstract 

The stability of static solutions of the spherically symmetric, 
asymptotically flat Einstein- Vlasov system is studied using a Hamilto- 
nian approach based on energy-Casimir functionals. The main result 
is a coercivity estimate for the quadratic part of the expansion of the 
natural energy-Casimir functional about an isotropic steady state. The 
estimate shows in a quantified way that this quadratic part is positive 
definite on a class of linearly dynamically accessible perturbations, pro- 
vided the particle distribution of the steady state is a strictly decreas- 
ing function of the particle energy and provided the steady state is not 
too relativistic. This should be an essential step in a fully non-linear 
stability analysis for the Einstein- Vlasov system. In the present paper 
it is exploited for obtaining a linearized stability result. 

1 Introduction 

The aim of the present paper is to investigate the stability of spherically 
symmetric steady states of the Einstein- Vlasov system against spherically 
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symmetric perturbations. The system describes, in the context of general 
relativity, the evolution of an ensemble of particles which interact only via 
gravity. Galaxies or globular clusters, where the stars play the role of the 
particles, can be modeled as such ensembles, since collisions among stars 
are sufficiently rare to be neglected. The particle distribution is given by a 
density function / on the tangent bundle TM of the spacetime manifold M. 
We assume that all particles have the same rest mass which is normalized 
to unity. Hence the particle distribution function is supported on the mass 
shell 

PM = {g a /3P a p = —(? and p a is future pointing} C TM. 

Here g a g denotes the Lorentz metric on the spacetime manifold M and p a 
denote the canonical momentum coordinates corresponding to a choice of 
local coordinates x a on M; Greek indices always run from to 3, and we 
have a specific reason for making the dependence on the speed of light c 
explicit. We assume that the coordinates are chosen such that 

ds 2 = c 2 goodt 2 + g a bdx a dx b 

where Latin indices run from 1 to 3 and t = x° should be thought of as a 
timelike coordinate. On the mass shell, p° can be expressed by the remaining 

coordinates, 

P° = a/^V 7 ! + c~ 2 g abP a p\ 

and / = f(t,x a ,p b ) > 0. The Einstein- Vlasov system now consists of the 
Einstein field equations 

G afi = 8vr C - 4 T^ (1.1) 
coupled to the Vlasov equation 

p°d t f + P a d x af - r^5„./ = o (i.2) 

via the following definition of the energy momentum tensor: 
rp | ,i/2 f t dp 1 dp 2 dp 3 

T a/ 3 = c\g\ 1 / PaPpf • (1-3) 

J —Po 

Here \g\ denotes the modulus of the determinant of the metric, and rjg are 
the Christoffel symbols induced by the metric. We note that the character- 
istic system of the Vlasov equation (jl.2p are the geodesic equations written 
as a first order system on the mass shell PM which is invariant under the 
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geodesic flow. For more background on the Einstein- Vlasov equation we 
refer to pQ. 

The system possesses a large collection of static, spherically symmetric 
solutions obtained for example in [181 E31 121] • The stability properties of 
these steady states have essentially not yet been investigated in the mathe- 
matics literature; numerical investigations of this question were reported on 
in [3], and in [25] the author employed variational methods for constructing 
steady states as minimizers of certain energy-Casimir type functionals. Ex- 
amples from the astrophysics literature where this stability issue is discussed 



In contrast, for the Vlasov-Poisson system, which arises from the 
Einstein- Vlasov system in the non-relativistic limit, considerable mathemat- 
ical progress on the question of the stability of spherically symmetric steady 
states has been made, cf. [T[ [HJ, O QUI CE31 E3 [H] and the references in the 
review articles [171 [20] . For isotropic steady states where the particle dis- 
tribution /o depends only on the local or particle energy the basic stability 
condition is that this dependence is strictly decreasing on the support of fo. 
In addition to the conserved energy or Hamiltonian the transport structure 
of the Vlasov equation gives rise to a continuum of conserved quantities, the 
so-called Casimir functionals. Using the latter the dynamics of the system 
can in a natural way be restricted to a leaf Sf of perturbations g : M 6 — > M. 
which have the same level sets as the steady state fo under investigation. On 
this leaf it is then possible to establish a coercivity estimate for the second 
variation of the Hamiltonian about fo- Such an estimate sometimes goes 
by the name of Antonov's stability bound. Under the assumption of spheri- 
cal symmetry this estimate was used in j9[ [15] to obtain non-linear stability 
against spherically symmetric perturbations for a large class of steady states 
/o- The case of general perturbations was finally completed in |16j . 




For the present case of the Einstein- Vlasov system we intend to employ 
an analogous approach. However, the numerical investigations in [H [27] 
indicate an essential difference in the stability behavior between the Vlasov- 
Poisson and the Einstein- Vlasov systems. For a given microscopic equation 
of state, i.e., a given dependence of fo on the particle energy, there typi- 
cally exists a one-parameter family of corresponding steady states. For the 
Vlasov-Poisson system all the members of one such family show the same 
stability behavior, but this is not so for the Einstein- Vlasov system. This re- 
markable phenomenon was first conjectured and numerically observed in the 
physics literature, most notably in the work of Ze'ldovitch [26\ 127], Here, 
a family of steady states with the same microscopic equation of state is 
parametrized by the central redshift which is a measure of how close to 
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Newtonian the steady state is. Provided that the microscopic equation of 
state is strictly decreasing, the steady states in such a family are stable 
when they are close to the Newtonian regime, i.e., when the central redshift 
is small, but they become unstable as this parameter increases beyond a 
certain threshold. For more precise statements of this behavior and the role 
of the so-called fractional binding energy in this context we refer to [1] and 
the original literature [26l [27] . 

Any stability analysis for the Einstein- Vlasov system must be able to 
reflect and deal with the above essential differences to the Newtonian case, 
and the present one does. We prove linear stability of a suitably defined 
one-parameter family of steady states with small central redshift, cf. The- 
orem 16.21 The main tool for the result is a quantified coercivity estimate 
for the second variation of the energy, analogous to the Antonov coercivity 
bound in the Newtonian case, cf. Theorem 14.21 This is a non-linear estimate 
that crucially depends both on the symplectic structure of the Einstein- 
Vlasov system and on properties of the Einstein field equations satisfied by 
the steady states. 

The paper proceeds as follows. In the next section we formulate the sys- 
tem in coordinates which are suitable for the stability analysis and explicitly 
introduce 7 = 1/c 2 as a parameter. We state the form of the ADM mass 
(or energy) and the Casimir functionals in these coordinates, and compute 
the quadratic term D 2/ Hc{fo) in the expansion about the steady state /o 
of a suitable energy- Casimir functional which is chosen such that the linear 
part in the expansion vanishes. In Section [3] we identify a class of linearly 
dynamically accessible perturbations which form the natural tangent space 
of the leaf Sf through /o which consists of all states which preserve all the 
Casimir constraints. The main result of our paper is shown in Section U) 
We prove that the second variation of the energy-Casimir is positive definite 
along the linearly dynamically accessible perturbations, provided that the 
parameter 7 = 1/c 2 is small enough. As in the case of the Vlasov- Poisson 
system this positivity result should play an important role in a future, fully 
non-linear stability analysis. In the present paper we restrict ourselves to 
drawing a conclusion on linearized stability. In order to derive such a result 
we analyze in Section [5] the linearized Einstein- Vlasov system. We provide a 
suitable existence and uniqueness theory for the corresponding initial value 
problem and show that the class of linearly dynamically accessible perturba- 
tions is invariant under the flow of the linearized system and that this flow 
preserves the quantity D 2 T-lc(fo)- In Section [6] we finally state and prove 
our result on linearized stability where we also deal with the fact that for 
our coercivity estimate the quantity 7 = 1/c 2 , which in a given set of units 
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is a specific constant, has to be chosen sufficiently small. This is done by 
exploiting a scaling symmetry of the problem which relates the Einstein- 
Vlasov system with unit speed of light c = 1 to the problem where the 
speed of light assumes a prescribed value c. In this analysis the smallness 
requirement for 7 = 1/c 2 translates into the requirement that the steady 
states under investigation are close to Newtonian. As explained above, such 
a sensitivity of the stability properties to being close to Newtonian or very 
relativistic is to be expected. 

Our analysis is restricted to spherically symmetric steady states /o and 
their stability against spherically symmetric perturbations. In particular 
the latter restriction is undesirable from a physics point of view. To remove 
it remains an important and highly challenging open problem; the existence 
of axially symmetric steady states was shown recently in [3]. 

To conclude this introduction we mention that the global existence result 
for the spherically symmetric Einstein- Vlasov system with small initial data 
[2T] can be considered as a stability result for the vacuum solution, but the 
techniques required for the stability analysis of non-trivial steady states are 
completely different from such small data results. Moreover, the question of 
weak cosmic censorship and black hole formation for the asymptotically flat 
Einstein- Vlasov system is studied in [6] and the latter question is addressed 
by different methods in [2J. A linearized stability analysis for the Vlasov- 
Poisson system is much easier than for the present case and was performed 
in [5] . Concerning the Hamiltonian approach for the Einstein- Vlasov system 
the reference |12j has been a most valuable source of inspiration. 

2 The spherically symmetric Einstein-Vlasov sys- 
tem and energy- Casimir functionals 

We consider a spherically symmetric and asymptotically flat spacetime. In 
Schwarzschild coordinates the metric takes the form 

ds * = _ c 2 e 2»(t,r) dt 2 + e 2X(t,r) dr 2 + ^2 + ^2 Q ^-g 

Here t £ R is the time coordinate, r £ [0, oo[ is the area radius, i.e., 4irr 2 
is the area of the orbit of the symmetry group SO(3) labeled by r, and the 
angles 9 £ [0, tt] and cp £ [0, 2ir] parametrize these orbits; c denotes the speed 
of light. The spacetime is required to be asymptotically flat with a regular 
center which corresponds to the boundary conditions 

lim \(t, r) = lim fj,(t, r) = = X(t, 0). 
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In order to formulate the Einstein- Vlasov system we write x = (x a ) = 
r (sin 6 cos (ft, sin 9 sin (p, cos 9) . Let p a denote the canonical momentum co- 
ordinates corresponding to the spacetime coordinates (x a ) = (t . X . X ■ 37 ^ 
and define 

v a =p a + (e A - 1)^ — , where x-p = 5 ab x a p b . 
r r 

This has the advantage that by the mass shell condition 

po = — e M -v/l + 7M 2 , where |i>| 2 = 5 a bV a v b and 7 := 
We introduce the abbreviation 

(v) := Vl + Tkl 2 - (2-2) 

In the canonical momentum variables the metric would also appear under 
the square root sign. The Einstein- Vlasov system now takes the following 
form: 

d t f + e^ A ^- ■ d x f - (A— + e^ x -p' (v)) --d v f = 0, (2.3) 

{v) \ r 7 / r 



e~ 2A (2rA'- 1) + 1 = 8irjr 2 p, 
e~ 2X {2rn' + 1) - 1 = 87r 7 2 r 2 p, 
-e~^~ A A = 47r7rj, 



e 

where 



-2A 



p(t,x) = J f(t,x,v) (v) dv, 
p(t,x) = J f(t,x,v)(^y ^, 

2 



q(t,x) = J f(t,x,' 



r 

xxv 



dv 

w 



(2.4) 
(2.5) 
(2.6) 



(V + (/*' " W + ;)) - (A + A(A - /i)) = Wq, (2-7) 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



Here d x f and d v f denote the gradients of / with respect to the x- and v 
variable respectively, and and ' denote partial derivatives with respect to 
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t or r respectively. Unless explicitly noted otherwise integrals extend over 
M 3 . Spherical symmetry means that 

f(t,x,v) = f(t,Ax,Av), x,v£R 3 , ,46 SO (3), 

i.e., / is invariant under the canonical action of SO (3) on the mass shell. 
As a consequence the spatial densities defined in (|2.8p - (|2,lip are actually 
functions of t and r = We refer to (TjJj for details on the derivation of 
this form of the equations. It has been used both in investigations of global 
existence for small data |21j and of the formation of black holes [2J. The 
former investigation also provides a local existence and uniqueness theorem 
for the spherically symmetric Einstein- Vlasov system for smooth, compactly 
supported and non-negative initial data 

/| t =o = /eC c 1 (]R 6 ), />0, (2.12) 

which are compatible with (|2.4p in the sense explained shortly. 

In a stability analysis the following aspect of the choice of coordinates 
should be noted: In order to compare unperturbed and perturbed quantities 
one must define some identification of points in the unperturbed with points 
in the perturbed spacetime where the various quantities are to be evaluated. 
A priori, no natural such identification exists, and for the Einstein- Vlasov 
system this question must be faced not only for points on spacetimes but for 
points on their tangent bundles or corresponding mass shells respectively. 
Throughout our analysis we identify points which have the same coordinates 
in the above set-up. Besides the fact that this choice works it is from a 
physics point of view motivated in [llj . 

For what follows it is important to note that given a spherically symmet- 
ric state / G Cc(K 6 ) we can explicitly and uniquely solve the field equation 
(|2,4p for A under the boundary condition A(0) = 0, 

e- 2 ^) = l- 7 ^, (2.13) 
r 

where 

rr 

2 



m{r) = 4-7T / s p(s) ds = / / (v) f(x, v) dv dx; 

JO J\x\<rJ 

we will occasionally write Xj, pj, rrif to emphasize that these quantities are 
determined by /. Clearly, for (I2.13P to define A on all of [0, oo[ we have to 
restrict the set of admissible states. We call a state / G C^(IR 6 ) admissible 
iff it is non-negative, spherically symmetric, and 

2m,f(r) 
7 — < 1, r > 0; 
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the initial data posed in (|2.12p must be admissible in this sense. 

We now define the ADM mass and the Casimir functionals as functionals 
on the set of admissible states: 



H(f) = J J (v)fdvdx, (2.14) 
C(f)= [[ e x f X (f)dvdx, (2.15) 



where x G C 1 (M) with x(0) = 0. These quantities are conserved along 
solutions launched by admissible initial data: 

n(f(t)) = n(f), c(f(t))=c(f) 

as long as the solution exists; for functions / = f(t,x,v) we write f{t) := 
f(t, •, •), and the analogous notational convention applies to functions of t 
and x or t and r. 

We need to recall some facts about steady states of the system (|2.3j) — 
(|2.1Up . a more detailed discussion including their dependence on 7 is given 
at the beginning of Section [H For p = po (r) time-independent and given 
we define the local or particle energy 

E = E(x,v) = e w,(r) (v) = e w(r yi + 7M 2 - (2-16) 

The ansatz 

f Q (x,v)) = <j>{E(x,v)) (2.17) 

with some suitable microscopic equation of state <j> = <ft(E) then satisfies 
the time-independent Vlasov equation and reduces the problem of finding 
stationary solutions to solving the two field equations (12. 4|) . (I2.5D where the 
source terms po and po are now functionals of fi = /io- A necessary condition 
for obtaining a steady state with compact support and finite ADM mass is 
that 4>{E) = for E > Eq where Eq > is some cut-off energy. We 
assume that <p £ C 2 (] — oo,Eq[) n C(M) is strictly decreasing on ] — 00, -E ] 
and such that a corresponding compactly supported steady state (fo, Ao,£to) 
with induced spatial densities po,Po and finite ADM mass exists. Examples 
of such functions are provided in [181 IZ3 [21] . 

We now discuss on a formal level the energy- Casimir approach towards 
stability for the Einstein- Vlasov system. We consider an energy-Casimir 
functional 

Hcif):=n(f)+C(f), (2.18) 



s 



and we expand about some given steady state /o of the above form: 

«o(/o + Sf) = Hc(fo) + DH c (fo)(5f) + D 2 n c (f )(Sf, Sf) + 0((<5/) 3 ). 

(2.19) 

A rather lengthy and non-trivial formal computation reveals the following. 
If the function \ which generates the Casimir functional (|2.15p is chosen 
such that 

X'(fo) = X'(HE)) = -E, i.e. x' = ~r\ 

then 

DHc(fo)(Sf) = 

and 

D 2 Uc(fo)(6f,6f) = l ff \^E)\^f) 2dvdx 

1 f'°° 

- — / e^°- A ° (2r/x' + 1) {5X) 2 dr. (2.20) 
27 Jo 

Here 5 A should be expressed in terms of 5f through the variation of f)2. 13j) . cf. 
(|3.3p below, and it should be observed that on the support of the steady state 
4> is strictly decreasing and in particular one-to-one. Since (j)'(E) = outside 
the compact support of the steady state /o, the first integral makes sense 
only for perturbations 5f which are supported in the support of the steady 
state. This is automatically the case for the class of linearly dynamically 
accessible states defined in Section [3l Just as in the case of the Vlasov- 
Poisson system the central difficulty in the stability analysis arises from 
the fact that the two terms in fj2.20f) are of opposite sign. For the Vlasov- 
Poisson system is has been shown in |13| that -D 2 % c (/o) is positive definite 
on a suitably defined class of linearly dynamically accessible states, and this 
fact has played a central role in the non-linear stability analysis in [9], [10] ■ In 
Section|4]we prove an analogous positivity result for the case of the Einstein- 
Vlasov system, and we believe that this should be a useful step towards a 
non-linear stability result. In this context the precise relation of D 2 n c (f ) 
to the energy- Casimir functional ()2. 18[) is important which is why we have 
gone into the issue of the expansion (I2.19p . However, in the present paper we 
only exploit D 2 H c (fo) and its properties to obtain a linear stability result. 
For this the relation (|2.19p is in principle irrelevant, the important issue 
being that D 2 n c (f ) is a conserved quantity along solutions of the linearized 
Einstein- Vlasov system, linearized about the steady state /o- This is shown 
is Section [5j The consequences of our results for non-linear stability are 
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under investigation, and in this context the expansion (|2.19p will have to be 
made mathematically precise. 

To conclude this section we recall the usual Poisson bracket 

{f,g}:=d x f-d v g-d v f-d x g (2.21) 

for two continuously differentiable functions / and gof i,i;£f 3 . We shall 
repeatedly use the product rule 

{f,9h} = {f,g}h + {f,h}g (2.22) 

and the integration by parts formula 

[f{f,g}dvdx = (2.23) 

which together with the product rule implies that 

{/, g}h dvdx = - \\ {/, h}g dv dx. (2.24) 



Since throughout, our functions are spherically symmetric it is sometimes 
convenient to use the following coordinates which are adapted to this sym- 
metry: 

II X - V 2 

r = \x\, w = , L = \x x v\ ; 

r 

a distribution function / = f(x,v) is spherically symmetric iff by abuse of 
notation, 

/ = f(r,w,L). 

It is also worthwhile to note that L is conserved along the characteristics 
of the Vlasov equation (|2.3p . i.e., due to spherical symmetry angular mo- 
mentum is conserved along particle trajectories. If / and g are spherically 
symmetric and written in terms of r,w, L, then 

{/> 9} = d r f d w g - d w f d r g. (2.25) 

Notice further that 

-{f,E} = e^^-- d x f - -e^ (v) - ■ d v f, 
7 {v) 7 r 



which should be compared with the Vlasov equation (|2.3p . 
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3 Dynamically accessible states 



Let (/o,Ao,/xo) with (|2.17p be a fixed stationary solution of the spherically 
symmetric Einstein- Vlasov system whose stability we want to investigate. 
We call an admissible state / non-linearly dynamically accessible from /q iff 
for all x G ^(R) with x (0) = 0, 

C(/)=C(/ ). (3.1) 

This property is preserved by the flow of the Einstein- Vlasov system. The 
aim of the present section is to develop a suitable concept of linearly dy- 
namically accessible states on which the second variation of He at /o is 
positive definite while the first one vanishes. Moreover, the set of these 
linearly dynamically accessible states should be invariant under the flow of 
the linearized Einstein- Vlasov system. Formally, taking the first variation 
in (13. ip . a suitable definition for bf to be linearly dynamically accessible 
should be that 

DC(f )(5f) = JJe Xo {x'(fo)Sf + x(fo)SX) dvdx = (3.2) 

for all x G C^E) with x(0) = 0, where 

5A = 7 e 2Ao — s 2 Sp(s)ds (3.3) 
r Jo 

and 

5p( r ) = 5p(x) = J (v) 5f(x, v) dv. (3.4) 

Condition (|3,2p is from an analysis point of view not practical to work and 
prove estimates with. In order to obtain a more explicit condition on 5f we 
first note a simple fact. 



Lemma 3.1 Under the assumptions made above, 



X{h)dv = - 1 e^ I x '(fo)<P'(E)^dv, 

(u) 



and hence 



DC(f )(5f) = JJe*°Mo) 



6f - 7 e^JA 



w 



dv dx. 
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Proof. We use the expressions E = e Mo (v) and L = \x x v\ 2 as new 
integration variables and observe that the integrand is even in v. Hence 

27re- w » E 

av = — — -= . clE dL 

^r 2 ^JE 2 - e 2 ^> (1 + 7 L/r 2 ) 

27re _W) 



d E \/E 2 - e 2 ^ (1 + 7 L/r 2 ) dL. 



An integration by parts now gives the result: 
x(/o) 



/ / , X '^(E)W(E)Je 2 - e^o (l + 1 -)dEdL 

JO JeW. XJ n \ V \ r / 



X mE))(/> (E) cfc 



-7e M / x'(/o)0'(^)^^. 



□ 

Using this result we find that a variation 5/ satisfies the condition (|3.2p . if 

2 

e Ao 5/ _ 7e Mo+Ao a = {h, /„} (3.5) 



for some spherically symmetric generating function h £ C 2 (R 6 ). This is due 
to the fact that 

' r x'(fo){h,f }dvdx = 



for any such generating function h. Since 5X appears in the second term, 
(|3.5p is still not suitable as a definition for 5f, but we have the following 
result: 

Proposition 3.2 Let 

w 2 r 

( 5/:=e -Ao {/i,/o} + 47r7 3 re 2w+Ao / (-S)T^-y <f/(E(x,v))h(x,v)wdv (3.6) 
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for some spherically symmetric generating function h G C 2 (IR 6 ). If the 
corresponding variation of A is defined by \3.3\) . then 

S\ = 4vrr7 2 e w+Ao / (j)'(E) h(x, v) w dv. (3.7) 



Hence 5f satisfies both \3. 5|) and \3.2i) . States of the form 113. 6\) are called 
linearly dynamically accessible from /o. 

For the proof the following auxiliary result is needed which will also be used 
elsewhere. 

Lemma 3.3 The following identity holds: 

g— Mo g— 2Ao-Mo 

cj)'(E)w 2 dv = (7p + Po) = -, — o — ( A o + A^'o) • 

7 47T7 z r ' 

Proof. We note that 

d v <j>{E) = 4>'{E)e^^- 

and hence 

-■d v 4>{E) = ct>'{E)e^. 
r \v) 

This implies that 

-Mo 

7 J r 



(t)'{E)w 2 dv = e — I ~ ■ d v 4>{E) (v) w dv 

(P(E) + (v)) dv 



7 J 

g— MO g— 2Ao-MO 

= — (7Po + Po) = ~ 47ry2r (A + Mo) 5 

the last equality is due to the field equations (|2.4p . (|2.5p for the steady state. 

□ 

Proof of Proposition 13.21 The variation of p induced by 5f takes the 
form 



5p = e- Xo J (v) {h, f }dv + 4vrr7 3 e Ao+2w J cf)'{E)hwdv j <f/(E) 
Hence Lemma 13.31 implies that 

Sp = e' Xo [ (v) {h, f }dv - 7 e w - Ao (A + p' ) [ <p'(E)hw dv 



w 2 dv. 



13 



We consider the first term on the right hand side and find that 
(v) {h,f }dv 

(v) (d x h ■ - d v h • (v) \ <j>'{E) dv 

7 e w [ d x h ■ v<p'(E) dv — // e M0 / d v h ■ -(1 + ^\v\ 2 )4>'{E) dv 



7e 



d x h ■ v<f)'(E) + 2^wh(j)'{E) + ef n' wh<f)" \E) (v) 



dv. 



Hence 



|x|<r 



2Ao ± 

2 



r 



+ e 



2A 7 



\x\<r 



d x h ■ vcf)'(E) + 2/j,' wh(j)'(E) 
+e^fj,' wh(j)"(E) (v) - (Aq + Mo) <P'(E)hw dv 
whcf)' '(E) dv 

{ l i' Q -\' Q )w4>'{E)-4>"(E)e^ l i' Q w(v) 



Q /-to— Ao 



(A + M ) 



/i <iu (fx, 



and since the term in brackets vanishes the claim is proven. □ 
Remark. For linearly dynamically accessible states as defined in Proposi- 
tion E21 

5f = 4>'(E) (e~ X °{h,E} +47T7 3 re 2M+A °^y J <j>'(E) h(x,v) w dv 
= 4>'(E) (e~ Xo {h, E} + 7 re^A W ~ 



(v)J ' 

in particular, Sf vanishes outside the support of /q and the integrals in 
(12^01) are well defined. 



4 The coercivity estimate 

The central step in our stability analysis is a coercivity estimate for the 
second variation of the energy. We recall the definition of this expression 
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which is also referred to as the free energy: 

A(5f):=D 2 HcUo)(Sf,5f) (4.1) 

1 f f p A ° 1 t'°° 

-{5ff dvdx- — e M ~ Ao (2rp' + l) (5X) 2 dr. 



2JJ \<P'(E)\" J ' 2 7 

The estimate will hold provided 7 is sufficiently small. We need to make 
precise the assumptions on the steady states under consideration and need 
to discuss their behavior for 7—7-0, i.e., their relation to those of the Vlasov- 
Poisson system 

dtf + v- d x f - d x U ■ d v f = 0, (4.2) 
AC/ = 4vrp, lim U(t,x) = Q, (4.3) 

\x\ — >OD 

p(t,x) = J f(t,x,v)dv. (4.4) 

Here U = U(t, x) is the gravitational potential of the ensemble, and p = 
p(t, x) is its spatial mass density. The Newtonian particle energy is given by 

E = E(x, v) = h,v\ 2 + U(x), (4.5) 
and under suitable assumptions on <3? the ansatz 

f (x,v) =$ (h v \ 2 + U(x)) (4.6) 



leads to spherically symmetric steady states. 
Assumption on $. Let $ € C(R) D C 2 (] - 00, 0[) with 

$ = 0on [0, 00 [ and < on ] - 00, 0[ 

be such that for a parameter v < the semilinear Poisson equation 

1 r 00 

^(r 2 U'(r))' = 4vr / §{E)J2(E-U(r))dE (4.7) 
r ' Ju(r) 

has a unique solution U E C 2 ([0,oo[) with the central value U(0) = v and 
the property that U(Rq) = for some radius i?o > 0. 

The ansatz (|4.6p necessarily leads to a spherically symmetric steady state 
where U = U(r), r = \x\, and reduces the static Vlasov-Poisson system to 
the equation (|4.7p : the corresponding steady state is supported in space in 
the ball of radius Rq about the origin. We also remark that the boundary 
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condition for the potential at spatial infinity has been replaced by the condi- 
tion that the potential vanishes at the boundary of the support of the steady 
state. This is technically advantageous, and the original boundary condi- 
tion can be restored by a suitable shift of the potential and a corresponding 
cut-off energy E$ < 0. A well known example where our assumptions hold 
are the polytropic steady states given by 

{ ' \ , E>0 

with k G]0, 7/2[; more examples can be found in [21]. 

For the Einstein- Vlasov system the particle energy is given by 

E = E{x,v) = e w,(r) (v) = e M(r Vl + 7M 2 - 
We fix a function $ as above and write = 7^0- By the ansatz 

f (x,v) = $ (^yi + 7M 2 e 7 "° - ^) = $ (~^) 

the static Einstein- Vlasov system is reduced to a single equation for uq, 
namely 

^o( r ) = ~, a* r 1 , , », ( irh^voir)) + — / s 2 gi {v {s)) ds J , 

1 ~ 1F7 Jo s 2 5 7 (^o(s))ds V r / 

(4-8) 

where <7~ and /i 7 are smooth functions determined by which are such that 
Po(r) = <? 7 (fo(r)) and po(r) = h^(vo(r)). The important point is that 

/•oo 

g 7 (u (r))^ / $(E)y/2(E-v (r))dE as 7 ^ 0. (4.9) 
./i'o(r) 

It should be noticed that the ansatz (|2.17p has been adapted to yield the 
proper Newtonian limit, and we have the relation </>(•) = <3?(- — ^). For the 
details of these arguments we have to refer to |23j ; here we collect only the 
information which we need for the stability analysis. 

Proposition 4.1 There exist constants 70 > and C > such that for 
< 7 < 70 the equation j^.ff| ) has a unique solution uq G C 2 ([0, 00 [) with 
fo(0) = v. The resulting steady state (/o, Ao, po = 7^0) satisfies the following 
estimates: ^ 

\x\ + \v I < C and v' Q > — on supp/o, 

and 

1 1 Po 1 1 00 > llPolloo, 1 1 Ao| |oo, 1 1 ^0 1 1 00 , ||^ol|oo<C. 
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Sketch of the proof. By the assumption on <3? the Newtonian potential 
U has a zero at some radius Rq, and since U is strictly increasing it is 
strictly positive for r > Rq. The structure of the right hand side of (|4.8p 
and in particular the limiting behavior (|4.9p implies that vq converges to U 
uniformly on bounded intervals as 7 — > 0, in particular, vq must also have a 
zero at a radius close to Rq for 7 small, and the rest follows; for details we 
refer to [23]. □ 

The content of the following theorem is the coercivity estimate for the free 
energy A on linearly dynamically accessible perturbations. 



Theorem 4.2 There exist constants C* > and 7* > such that for any 
< 7 < 7* and any spherically symmetric function h £ C 2 (IR 6 ) which is 
odd in the v-variable the estimate 



A(5f)>C J J W(E)\ ((rwf 



E,— 

rw 



+ 7 2 |/i| 2 dvdx 



holds. Here 5f is the dynamically accessible perturbation generated by h 
according to \3. 6\) . 

Before proving this theorem we collect some auxiliary results. 



Lemma 4.3 Let h G C 2 (R 6 ) be spherically symmetric. Then the following 
estimate holds: 



2 -2A -Mo 

b'(E)\whdv) < (A'o + ^'o) / \4>'{E)\h 2 dv. 



Proof. By the Cauchy-Schwarz inequality, 



\(f)' (E)\wh dv ) < / \(j)'{E)\w 2 dv / \(j>' {E)\h 2 dv 



and applying Lemma 13.31 completes the proof. 



□ 



Lemma 4.4 The following identities hold: 



{E, rw} = e^rfiQ (v) - e Mo (v) + 



W 



{E,{E,rw}} 



rf4 + Mo + 7% 
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Proof. The first identity follows easily from the product rule (|2,22p and the 
definition (|2.16p of E, for the second identity we use the first one and again 
the product rule. □ 
Proof of Theorem 14.21 Before starting our estimates we mention that 
throughout the following proof, < 7 < 70 as in Proposition 14.11 and C 
denotes a positive constant which may change its value from line to line, 
depends on the bounds provided in Proposition 14.11 and is independent of 
7- 

Step 1: Splitting A. Using the formula (j3.6|) for 5f and the definition (|4.ip 
of A, we expand the square of Sf appearing under the integral sign and 
rewrite A(5f, S\) in the following form: 

2A(5f)=A 1 (5f)+A 2 (5f), (4.10) 

where 

M6f) := // ^ \i E > h }\ 2dvdx 
1 f'°° 

__/ e^- Xo (2rfi' + l)(5X) 2 dr, 
7 Jo 

=:A u (5f) + A 12 (5f), (4.11) 

2 



A 2 (5f) := -2 7 J I \fi(E)\{E ) h}5\e> 10 ^dvdx 

+7 2 // |0'(£)|e 2 ^ +Ao ^((U) 2 (Wx 
J J (v) 

=:A 2 i(Sf)+A 2 2(Sf). (4.12) 

The idea behind this splitting is that .4,1 will yield the desired lower bound 
while A 2 is of higher order in 7 and can eventually be controlled by the 
positive contribution from A\. 

Step 2: The term A\ . Let us define 

n := — h. 

rw 

The function rj is well-defined for w = since h is odd in v. By the product 
rule (1232]) . 

{E,h} = rw{E,rj} +r){E,rw}, (4.13) 

and by QHgD , 

\{E, h}\ 2 = {rw) 2 \{E, r]}\ 2 + {E, r] 2 rw{E, rw}} - r] 2 rw{E, {E, rw}}. 
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Inserting this into An, we arrive at 

~ Xo \0'(E)\\{E,h}\ 2 dvdx = J J e- Xo \ft(E)\(rw) 2 \{E,T]}\ 2 dvdx 
+ [[ e- Xo \<j)'(E)\{E,r] 2 rw{E,rw}}dvdx 
e~ x °\<f)'(E)\ r] 2 rw{E, {E, rw}} dv dx 



=:U + V + W. (4.14) 

The term U is a positive definite contribution to A\. The terms V and W 
are more delicate since they will have to be combined with the negative term 
A\2- For the term V we use the formula (|2.24p and integrate by parts in v 
to arrive at 



V = — J J e~ A °{/o, t] 2 rw{E, rw}} dv dx 
= J J fo{e~ Xo ,r] 2 rw{E,rw}} dv dx 
= J J fod x (e~ x °) ■ d v (r] 2 rw{E, rw}}) dv dx 

d v fo • ^ e" A ° (-X' )r] 2 rw{E, rw} dv dx 

b'(E)e^^- ■ -e- Xo X' Q r] 2 rw{E,rw} dv dx. 

{v) r 

Using the first identity in Lemma 14.4} 

V = ~ 7 II \^( E )\ e2 ^ X ° X o^f (Vo («> " («> + A") dvdx 
= -jJJ W{E)\e 2 ^h 2 (V -^ + 
By the second identity in Lemma 14.41 

W = 1 J J \<j)\E)\e 2 ^- x »h 2 ^ + ^ + dvdx, 



— iT ) dv dx. 



r (v) 2 



and hence 



V + W = 1 JJ ^'(i^o-Ao^ _ x y Q + f^±K + ^L*t>) cfo dx. 



r (v)' 

(4.15) 
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In order to bound the expression A12 from below, we use the formula (|3.7p 
for SX and Lemma |4,3I 

7 Jo 



1 

7 Jo 

> -47T7 3 



(2r//o + l)16vr 2 r 2 7 4 e 2Ato+2Ao / <j>' \E)hw dv dr 



= 3/io+Ao 



(2^0 + 1)- 



— 2Aq— ho 
47T7 2 r 



(A'o + Mo) / |0'(^)|/i 2 ^dx 



- 7 / / ^(JEOIe^-^fc 2 ( 2 Mo A + 2(^) 2 + <7r ,/,-. ! I.I G) 



i,From (|4.15p and (|4.16p we have 



V + W - - 

7 Jo 



IMi— Aq 



(2r^ + l)(5A) 2 c?r 



> 7 // |0 / (E)|e 2 ' io - Ao ^ 2 



- 3// A' - 2( M ' ) 2 + 



(4.17) 



In order to estimate the term in the rectangular brackets on the right-hand 
side of (|4.17p . we express fi' Q ' via the field equation (|2,7p for the steady state 
and obtain, using the fact that the source term go in the latter equation is 
non- negative, 



[. . . J = 4tt7 q e - 2// A - 3(/z ) - — + — + 2— -5 — ^ 

r r r (u) r (v) 



> -2»>M + A ) - (/i ) 2 + ^0 



1 1 + — ( 1 

r 



1 



o 2 2A // , \ 1, l\1 , 7^0 1-7 ^ , A 7M 

= -8vr 7 re °i/ (po + 7Po) - 7 K) + 771— + ~ "77: 

r («) r ( w ) 

The observation that 

X'o = A' + // - // = 7 (4^r(p + 7Po)e 2A ° - ^) > -7^0 
implies that 

V ° 7 2 ( ^ru'^{p + 7 Po) + K) 2 + 2^ b|2 



■■■ >7 



r («)' 



7^72 [l-7 f8vrr 2 (w) 2 e 2A °(p + 7Po ) + r («) 2 ^ + 2| V | 2 ) 
r {v) L v ' 
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Using Proposition 14.11 there exists a constant 71 e]0, 70] such that for all 

< 7 < 71, 



[•••]> 



7^o 
4 r ' 



(4.18) 



Therefore, from (|4.14p and (|4.16p . assuming this smallness condition for 7 
we obtain the crucial estimate 

Ai(df) > J J \(j)'{E)\ (e~ x °(rw) 2 \{E,r)}\ 2 + ^e^-^^h 2 ) dv dx 

>C if \^'(E)\((rw) 2 \{E,r]}\ 2 +j 2 h 2 )dvdx, (4.19) 



where we again used the bounds from Proposition 14.11 

Step 3: The term A2- Since this term is of higher order in 7, we expect 
to be able to prove that it is small compared to A\ for suitably small 7. 
Using the decomposition (|4,13p and keeping in mind that 77 = h/rw and the 
formula (|3,7p for 5\ we can rewrite the part A21 as follows: 

A21 = ff \(t>'{E)\e 2 ^ 0+Xo ^-{E,h} (J (j)'{E)hwd%\dvdx 

= -8tt~/ 3 JJ \(f)\E)\e 2 » 0+Xo ^-{E,r)} (J (j)'{E)hw dv^j dv dx 
JJ \4>'(E)\e 2 ^ +Xo ^{E,rw}h (J ^'(E)hwdv^j dv dx 

=: X + Y. 
By Lemma 14.31 



-Mo 



7 



C 



7 



\<P'(E)\whdv) < — — (p + 7Po) / \4>'{E)\h 2 dv < - I \4>'{E)\h 2 dv. 

(4.20) 

In addition, Lemma f3.3l together with the bounds from Proposition ^. 1 l imply 
that 



sup / \(j)'{E)\w 2 dv < 



C 
7 ' 



(4.21) 



The bounds from Proposition 14.11 the Cauchy-Schwarz inequality, and the 
estimates P~20]) . (jOTj) imply that 



\X\ < C7 3 



\4>'(E) 1 1/2 w\ 4f (E) 1 1/2 rw{E, 77} dv 



\<t>'{E)\hw dv 



dx 



< C 7 5/2 ( / / \(j>'(E)\\rw{E,r]}\ 2 dvdx 



1/2 



\(j>'(E)\h 2 dvdx 



1/2 
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Hence by (ETT9|) . 

\x\ < cvii. 

In order to estimate the term Y we rewrite the first identity of Lemma 14.41 
{E, rw} = e^nu'r, (v) + e Mo f- (v) + -^-) = je^ (rv'n-^x 

V («>/ V w 

Thus, 

sup |{£", ru;}| < C7. 

supp /o 

Using this together with the estimates from Proposition 14.11 we proceed as 
above to find that 

\Y\ < C~/ 7/2 (^jj \<J)'(E)\h 2 dvdx^ (J- J J \cj>'{E)\h 2 dvdx\ ' 
< C~iA x . 

/,From the above estimates for \X\ and \Y\ it follows that 

A 2 > -\X\ - \Y\ > -C7A1. 
^From this estimate and (|4.10p . we finally infer that 

1 1 

A > -Ai - CjAi > -Ai 

provided 7 is sufficiently small. In view of (|4.19p the proof is complete. □ 

Remark. The estimate (|4.19j) shows that an equivalent way of stating our 
coercivity estimate is 

A(Sf) >\JJ \^'{E)\ (e- x »{rwf\{E, n }\ 2 + f^Mo-Ao^^ dy ^ 

which should be compared with the coercivity estimate in the Newtonian 
case, cf. {9} Lemma 1.1]. 

If we want to use Theorem l4.2l to deduce a stability result the restriction 
that h is odd in v is a problem because the generating functions of our 
perturbations need not be odd and more importantly, even if they were this 
property is not preserved under the linearized flow. However, the restriction 
is easily removed. To this end, we define for a function h = h(x, v) its even 
and odd parts with respect to v as usual by 

1 1 

h + (x,v) := - (h(x,v) + h(x, —v)) , h^(x,v) := - (h(x,v) - h(x, —v)) . 
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Corollary 4.5 For any spherically symmetric function h 6 C 2 
timate 



the es- 



A(Sf)>C* J J W{E)\ Urw) 



h. 



rw 



dv dx 



+ 



" Ao \</>'(E)\ \{E,h + }\ 2 dvdx 



holds. Here 5f is the dynamically accessible perturbation generated by h, 
C* > and 7* > are as in Theorem \4-%\ and < 7 < 7* . 



Proof. We split h into its even and odd parts, h = h+ + h-. Then 

5\ = 47rr7 2 e M0+A ° f <j>'(E) h-(x,v) w dv, 



since h + does not contribute to the last integral. Since /o is even in v this 
implies that 

Ao j 



6f. = e~^{h + Jo}, 



w 



Sf+ = e- Xo {h^J } + 7 e^ ( / ) '(E)—5X, 

in particular, the even part of 5f is the dynamically accessible perturbation 
induced by the odd part of h. Hence 



A(5f) = A(Sf+) + J J 



6 wm 2 



" A o I J/ 



|{£,/i + }| 2 (focte, 



since the integrand of the mixed term is odd in v. The assertion follows if 
we now apply Theorem 14.21 □ 



5 The linearized Einstein- Vlasov system 

In order to linearize the spherically symmetric Einstein- Vlasov system (12 . 3ft - 
(|2.10p about a given steady state (/o, Ao,^o) we write 

/(*) = fo + Sf(t), A(t) = A + SX(t), n(t) = fio + Sfi(t), 

substitute this into the system, use the fact that (/o,Ao,A*o) is a solution, 
and drop all terms beyond the linear ones in (Sf,SX,Sfj,). At the moment it 
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may not be obvious that this notation is consistent with the one from the 
previous sections, because the 5 terms now have a different meaning. But 
it turns out below that the notation is indeed consistent, and we obtain the 
following system: 

d t 5f + -e~ x °{6f,E}- (~f5\w + e>* o - Xo 5fi'(v))0'(E)—— = O, (5.1) 

7 V /(f) 



e~ 2X ° (r5X' - SX {2rX' - 1)) = 4irjr 2 dp, (5.2) 

e- 2Xo {rSp! - 5X (2r// + 1)) = 4vr 7 2 r 2 5p, (5.3) 

where 

Sp(t,r) = J (v)Sf(t,x,v)dv, (5.4) 

5p(t,r) = J —5f(t,x,v)dv. (5.5) 

Since the condition of a regular center implies that 5X(t, 0) = the quantity 
5X is determined by (|5.2|) as 

SX(t,r) = 7 e 2Ao — / s 2 5p(t,s)ds, (5.6) 



which agrees with the previous definition of 5X in (|3.3p . The question arises 
whether the linearized version of (|2.6[) follows from the linearized system 
stated above; we will make no use of the linearized version of (|2.7p to which 
this question of course applies as well. Indeed, 

SX = -4^7re w+Ao (5j, (5.7) 

where 



$j{t, r ) = J w 5f(t,x,v) dv. (5i 
To see this we observe first that by (|5.6p . 



SX(t, r) = 7 e 2Ao — / s 2 d t 6 f >{ l.s) r/.*. 



By (15. ip and integration by parts, 

d t 5p = -e^- Xo j v ■ d x 5f dv - 2p' e^- x ° j wSfdv- SX (p + 7Po)- 
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If we substitute this into the equation for 5 A and integrate the first term by 
parts with respect to x we find that 



r5\ + AiTr 2 ~{e^ +Xo 5j = -4irje 2X ° 




Since the term in parenthesis vanishes at the origin it vanishes everywhere, 
which is the assertion. 

Since we are restricting our linear stability theory to the class of dynami- 
cally accessible perturbations, we limit ourselves to an existence theorem for 
such solutions. In particular, we show the essential fact that linearly dynam- 
ically accessible data retain this property under the flow of the linearized 
system. 

In order to avoid purely technical complications we assume from now on 
that in addition to the previous assumptions, $ £ C 2 (M). We comment on 
this assumption in a remark at the end of this section. 

Theorem 5.1 Let h £ C 2 (R 6 ) be a spherically symmetric function which 
according to \3. 6\) generates a linearly dynamically accessible perturbation 
5f. Then there exists a unique solution 5f £ C 1 ([0, oo[x]R 6 ) to the linearized 

o 

Einstein- Vlasov system 115. l\) - [5l)\) with (5/(0) = 5f. Furthermore, there 
exists h £ C 1,2 ([0, oo[xM ) such that 



Sf(t) = e- Xo {h(t),f Q } + 4ir- f 3 re 2flo+Xo <i)'(E) T - / ^'(E)hwdv, (5.9) 



i.e., 5f(t) is linearly dynamically accessible. The generating function h is 
the unique solution to the transport equation 



with initial value h(0) = h 

Proof. The proof proceeds as follows. First we establish existence, unique- 
ness, and regularity of the solution h to the equation (|5.10p . where 5X and 
6 fi are defined via the field equations (|5.2p and (I5.3p . with source terms 
induced by 5f as defined in (|5.9p . Then we prove that (5f,5\,S/j,) indeed 
solves the linearized Einstein- Vlasov system. 

Equation ()5.10j) is a first order inhomogeneous transport equation and 
its characteristics s h-» (X(s, t, x, v), V(s, t, x, v)) are defined as the solutions 





(5.10) 



o 
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of 

~ _ Hq(x)-\q(x) V 

(v) ' 

y = _I e MW-AoW V/io(x ) ^ 

7 

with initial condition 

X(t, t, x, v) = x, V(t, t,x,v) = v. 

The metric coefficients of the steady state can of course be viewed as func- 
tions on R 3 , and as such, A G C 2 (R 3 ) and /x G C 3 (M 3 ), cf. [231 [M]- Using 
the abbreviations z = (a:, v) and Z = (X, V) we see that Z(s, t, •) : M 6 — > M 6 
is a C 2 diffeomorphism. Now assume that h is a solution to (I5.10p . Then 
integration along the characteristics implies that 

h(t,z) = h(Z(0,t,z))- [ (ef^SX^-r + -E5u) (s,Z(s,t,z))ds. (5.11) 

Jo V W 7 / 

We wish to find a solution to this equation, where 5X and <5// are given as the 
solutions to the field equations (|5.2I) and (|5.3f) with source terms induced by 
Sf which in turn is defined via (|5.9p . 

To construct solutions to (|5.1ip . we apply a simple iteration scheme. For 
any spherically symmetric function g G C 2 (R 6 ) we define the function 5f g 
by (|5.9p with h replaced by g. Clearly, Sf g G C 1 (IR 6 ), and Sf g is supported 
in the support of fo- The induced source terms Sp g and 5p g have the same 
regularity and are compactly supported as well. The equation (15. 6p shows 



that the induced metric component 5X g G C 2 (M 3 ). Moreover, the formula 
in Proposition 13.21 shows that 5X g is also compactly supported. Next we can 
define 8n' g by (|5.3p and using the boundary condition <5// s (oo) = we find 
that 5[i„ G C 2 (M 3 ) is compactly supported as well. Moreover, 



5\ g \\c2 + \\^ g \\ C 2 < C [\\5 Pg \U + \\5 Pg \\ ci < C\\5f g \\ ci < C\\g 



\c. 



(5.12) 

where the constant depends on the given steady state and the norms extend 
only over the support of the steady state which is compact. Suppose now 
that g G C 1 ' 2 ([0, oo[xM 6 ) so that g(t) G C 2 (M 6 ) for t > 0. Motivated by 
(|5.1ip we define 

(Tg)(t,z) :=&(Z(0,t,z))- J* ^ ^ 9 (s)^ + ^.)) (Z(a,t,z))da. 

(5.13) 



26 



It is straight forward to see that Tg G C 1,2 ([0, oo[xlR 6 ), and the linearity of 
the problem together with (|5.12[) imply that for 51,52 £ C 1,2 ([0, oof, M 6 ) we 
obtain the estimate 

\\T 9l (t)-Tg 2 (t)\\ c? <C f t \\g 1 (s)-g 2 (s)\\ C 2. (5.14) 







Now let 

o 

h (t,z) := h(z), h n+ i :=Th n . 

The estimate (|5.14|) implies that h n converges i-locally uniformly to some 
h € C 1,2 ([0, oo[xM 6 ). Using again the estimates in (|5,12p and linearity it 
follows that Sfh n converges t-locally uniformly to 5fh E C 1 ([0, oo[xM 6 ) and 
5Xh n , 5fih n converge t-locally uniformly to 5Xh,5/j.h G C 1,2 ([0, oo[xlR 6 ). In 
particular, h solves (I5.11j) and hence also (15.10p . The uniqueness of the 
solution is clear. 

In order to show that (5f, SX, Sfi) solves the linearized Einstein- Vlasov 
system, it only remains to check the Vlasov equation (|5.ip as the two field 
equations (|5,2p and (|5.3p hold by definition of 5 X and 5fi. The definition of 
5f in terms of h implies that 

2 

d t Sf = e~ x °{d t h, f } + ^SXtfiE)^. 



Comparing this with (|5.ip we see that the latter equation becomes equivalent 
to the relation 

{d t h, /„} = --{5f, E} - e 2 ^5^'w <j>'{E) 
7 

= _I {e -Ao {/i;/o})E} _ \e^5X<P'(E)^,E) -e 2 ^5^w4>'(E). 
7 I W J 

The fact that / = 4>(E) and {(j)'(E), E} = = {E, E} together with the 
product rule (|2.22p imply that this relation is again equivalent to 

{dthJo} = --{e- x y'(E){h,E},E} - (e^SX 4>'{E)% e) 
7 I (v) J 

-e^dfi'wftiE) 

= ~cj>'(E){e- x °{h,E},E} - 4>'{E) L^5X^e\ 
7 I (v) J 

--4>'(E){E5fi,E}, 

7 
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i.e. 



d t h + -e- Xo {h, E} + e^SX^- + -ESp, fo } = 



7 ' (v) 7 

which is implied by (|5.10p . □ 

Remark. We emphasize that a global existence and uniqueness result for 
general smooth data holds as well. The proof follows a simple iteration 
scheme analogous to the corresponding result for the Vlasov-Poisson case [5] . 
We also note that the values of h outside the support of fo are actually 
irrelevant and it would be sufficient to consider the generating function h as 
defined only on this support. 

The free energy A defined in (|4.ip is conserved along solutions of the 
linearized system. This fact, which is clearly important for our stability 
result, is shown next. 

Proposition 5.2 Any linearly dynamically accessible solution as construct- 
ed in Theorem \5. 1\ preserves the energy A. 

Remark. The above assertion is true for any sufficiently smooth solution 
of the linearized system, provided in particular that the first integral in 
A(5f(t)) is defined. As the remark at the end of Section [3] shows, this 
is the case for linearly dynamically accessible solutions as constructed in 
Theorem 15.11 In the following proof we make use of this structure only to 
guarantee the existence of the otherwise questionable integrals. 

Proof of Proposition 15.21 Clearly, 
I A ff e *ojf dvdx= ff e \°V_ dtSfdvdx 



2dtJJ W(E)\ J J W{E)\ 



{6f, E} dv dx 



1<P'{E) 

an , 

5f (jSXw + e^-^Sfi' {v} ) dv dx 



(V) 

5f 



{8f, E} dv dx 



/(E) 

je Xo+flo 5p5X + e 2fM0 5j5fi') dx. 



For the first term on the right hand side we use the product rule (|2.22p and 
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the identity (|2.24p to conclude that 



{E,5f}dvdx = -- // - jr -r{E i 5r}dvdx 



V(E)\ ' JJ 277 
1 



2 77 ^( B '^*" b - a 



Hence the linearized field equations (15. 3D and (15. 7h imply that 

Id /"/"-An */ 

oMo— Ao pMo-A 



/ / e Xo -^—7dvdx = - f (je Xo+,10 5p5X + ^SjStA 
2 <ii 77 |0 '(E) | 7 v ' 

= - 1 ( 6X — (»V ~ 6X ( 2 ^o + !)) - — 6X ) r " dr 



-I /-oo 

= _ / e Mo-A (2r//(, + 1) 6X S\ dr 
7 7o 

= 2^Jt J eMO " A °( 2 ^o + 1 )(^) 2 dr 

as required, and the proof is complete. □ 

Remark. In Theorem I5.1| the assumption <I> 6 C 2 (IR) can be relaxed to 
a C 1 assumption. The existence theory can then be developed for h € 
C 1 ([0, oo[xR 6 ). Since 5f is then only continuous, it becomes more technical 
to justify Proposition 15.21 but this can be done analogously to the proof of 
[5j Theorem 5.1]. 



6 Linear stability 

Our coercivity estimate in the form of Corollary 14.51 and the conservation of 
the free energy according to Proposition 15 . 21 immediately imply the following 
result. 

Theorem 6.1 Let satisfy the assumptions stated above, let 7* and C* be 
as in Theorem and let < 7 < 7* . Then the corresponding steady state 
introduced in Proposition \4-l\ is linearly stable in the following sense: For 
any spherically symmetric function hE C 2 (M 6 ) the solution of the linearized 

_ i o 

Einstein- Vlasov system S5. ij) - l5T5j) with the dynamically accessible state 5f 

o _ 

generated by h according to S3. 6\) as initial datum satisfies for all times t > 
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the estimate 



C" // \<p'(E)\ \(rw) 



rw 



+ 7 2 |/i_(t)| 2 dvdx 



e- Xo \<P'(E)\ \{E,h + (t)}\ 2 dvdx<A(5f). 

Remark. The left hand side is an acceptable measure of the size of the 
linearized perturbation from the given steady state, in particular, if it van- 
ishes then so does the perturbation 5f(t). The fact that the left hand side 
controls h{t) only on the support of the steady state is natural in a linearized 
approach. The right hand side can clearly be made as small as desired by 
making the initial perturbation small in a suitable sense. 

The result above is acceptable as a linear stability result, except for the 
fact that 7 = 1/c 2 has to be chosen small when in a given set of units this 
quantity has a definite value. We therefore recast our result into one on the 
stability of a one-parameter family of steady states of the Einstein- Vlasov 
system with 7 = 1 by using the scaling properties of the system. This will 
also provide a quantitative version of the Ze'ldovitch observation which was 
explained in the introduction. 

In order to do so we now denote the Einstein- Vlasov system (I2.3p -( J27TT]1 . 
i.e., including the parameter 7, by (EV 7 ) and by (EV) the system with 7 = 1. 
A straight forward computation shows that if / is a solution of (EV 7 ) with 
the corresponding metric coefficients A and /i, then 

Tif(t, x, v) := 7 " 3/2 /(i, 7~ V2 x, l^v) 

defines a solution of (EV) with metric coefficients 

Ty\(t,r) := A(t j7 - 1/2 r), T>(t,r) := ^t^r). 

For a fixed (ft as specified above, the steady state (/o,Ao,7o) of (EV 7 ) with 
< 7 < 70 — this is actually a family of steady states, one for each (EV 7 ) — is 
now mapped into the one-parameter family of steady states of (EV) given 
by 

]0,7o] 3 7^ (fl AM = (TVo,^A 0! rVo); 

here 70 is from Proposition 14.11 It should be carefully noted that all the 
members of this family are steady states of (EV), i.e., of the system with 
7 = 1, and 7 G]0, 70] is now the parameter which parametrizes this family. 
In particular 

/o 7 = 1~ 3/2 HE) = ME) 
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where 

E(x,v) := e^ f Vl+W 2 and ^ := 7~ 3/2 </>. 

In order to understand what it means to move to smaller values of 7 in this 
family we observe that by construction fi^(W) = holds for a unique radius 
W > which determines the boundary of the spatial support of the steady 
state, i.e., the boundary or surface of the galaxy or globular cluster. At the 
center r = we have A*q (0) = ju The expression 

e ^ ( R1 ) 1 
Z := yjt: 1 = — 5 - 1 > 

is the redshift of a photon which is emitted at the center and received at 
the surface of the mass distribution, and it is a measure for how strong 
relativistic effects in the configuration are. Expressed in units where the 
speed of light c = 1 the above limit 7—7-0 means that we consider steady 
states for which z is close to 0, i.e., relativistic effects are weak. 

In order to recast our stability result for the case of (EV) we have to 
check how the various quantities which are involved behave under the scaling 
operator T 7 . To this end we rename the particle energy in the context of 
(EV 7 ) as 

£ 7 = E y (x,v) = e^Vl + TM 2 = £(7 1/2 £,7 1/2 v). 

In order to understand the behavior of the Poisson bracket (|2.2ip under the 
scaling operator T 7 , we define for a given function h = h(x,v), 

h^{x,v) := 7 - 1 / l (7 1 / 2 x,7 1/2 i»). 

Then the relation 

{hJJ}(x,v)=T^{h\f } 
holds. Next, a simple change of variables argument implies that 

Airre^ +X *4>'JE)— / <j)'(E)hwdv 

= T 7 ^47r 7 3 re w+A( y(£ 7 )^ J <f>'(E 7 )h^wds \ . 



In particular, if 



2 [■ 

'o{/ i ,/ 7 } + 4vrre' 1 o+^^( jB )^_ J <fi^(E) hwdv 
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is the linearly dynamically accessible state for (EV) generated by h, then 



Sf h = T^8f h ,, (6.2) 

where 5fhi is the associated linearly dynamically accessible state for (EV 7 ), 
generated by K 1 as defined in (|3.6|) . 

Finally, to understand the behavior of the the free energy we first 
must again be careful with the notation. The free energy associated with 
the steady state (/o,Ao,/io) of (EV 7 ) and defined in (|4.ip is still denoted 
by A, while the corresponding quantity associated with the steady state 
(flXluZ) of (EV) is defined by 

Ai{5f) := \ is \^M {5f)2dv dx ~ir e ^~ A ° ( 2r ^ + ^ {5x)2 dr ~ 

A simple calculation shows that 

A(5f) = j- 3 / 2 A~<(T^5f). (6.3) 

We now recast Theorem 16. 1\ which is a result for (EV 7 ), into a result for 
(EV). 

Theorem 6.2 Let <I> satisfy the assumptions stated above and let 7* and C* 
be as in Theorem \4-2\ Then provided < 7 < 7* the corresponding steady 
state (/q ', Aq, {j,q) of (EV) is linearly stable in the following sense: For any 
spherically symmetric function h£ C 2 (M 6 ) the solution fh of the linearized 

o 

Einstein- Vlasov system with the dynamically accessible state 5f generated by 
h according to \6.1\) as initial datum satisfies for all times t > the estimate 



c // WJE)\ H 



rw 



|/i_(t)| 2 dvdx 



+ \ jje- X o\^(E)\ \{E,h + (t)}\ 2 dvdx<A\5f). 

Proof. In the statement of the theorem the linearized Einstein- Vlasov sys- 
tem is now the system ()5.1|) — (|5.5|) with 7 = 1, and A' 7 is a conserved quantity, 
in particular along solutions which are linearly dynamically accessible from 
the steady state (/J ', Aq, /^o) about which the system (EV) was linearized. 
Hence by E2D and (IQl) . 



AH8°f) = A*mV)) = A\T^5f w {t)) = 7 3/2 'Ai5fh.it)). 
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If we apply Corollary 14.51 it follows that 



A^(5°f) > C* 7 3/2 



+ 7' 



.3/2 




Now we apply a change of variables to turn /i 7 into h under these integrals 
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